Impact of network structure on a model of diffusion and competitive interaction 
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We consider a model in which agents of different species move over a complex network, are subject to 
reproduction and compete for resources. The complementary roles of competition and diffusion produce a 
variety of fixed points, whose stability depends on the structure of the underlying complex network. The survival 
and death of species is influenced by the network degree distribution, clustering, degree-degree correlations and 
community structures. We found that the invasion of all the nodes by just one species is possible only in 
Erdos-Renyi and regular graphs, while networks with scale-free degree distribution, as those observed in real 
social, biological and technological systems, guarantee the co-existence of different species and therefore help 
enhancing species diversity. 
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In the last few years, complex networks have been the sub- 
ject of an increasingly large interest in the scientific commu- 
nity JlJ-llll- This is due to: i) the wide variety of complex 
systems that can be described as graphs with a complex topol- 
ogy; ii) the recent observation that various dynamical pro- 
cesses taking place on networks, such as epidemics JUl, ran- 
dom walks llJfl], synchronisation J3l and self-organised crit- 
ically HI, can be affected by the underlying network struc- 
ture. Concerning social interactions, many models have been 
proposed in the last decades to study evolution of relation- 
ships, culture segregation, opinion formation and propagation 
of new ideas by means of majority rules, melting or mi xing , 
and have also been studied on complex topologies Ill0l4l4ll . 
Particularly important contributions to the understanding of 
social dynamics taking place on networks have been provided 
by evolutionary game theory lfl5l [Trill . Usually, when an evo- 
lutionary game is studied on a graph, each individual is asso- 
ciated to a node of the graph, and the social relationships are 
represented by the links. Individuals interact with their neigh- 
bours on the graph by playing various evolutionary games, 
and different collective behaviours emerge, such as global 
cooperation or selfishness, accordin g to the structure of the 
underlying network of relationships 11171 [Till . These models 
on complex networks fail to catch one of the most important 
characteristic of real evolutionary systems, namely the pos- 
sibility for the individuals to move through a complex envi- 
ronment and to interact with a neighbourhood which changes 
over time lfT9tl . Some recent works have proposed extensions 
of evolutionary game theory to moving agents J20I |2lll . In 
these latter models, the individuals move over a continuous 
space or a discrete lattice, and play games with other individ- 
uals in their spatial neighbourhood. However, the hypothesis 
of a homogeneous and continuous space is too simplistic, and 
does not correspond to the structure of real social and techno- 
logical networks fTl-[3i l22[l . Metapopulation models with het- 
erogeneous connectivity patterns, which incorporate mobility 
over the nodes, local interaction at the nodes, and a complex 
network structure, have been recently proposed only in spe- 
cific contexts such as epidemic spreading f23ll or chemical 
reactions l24ll . In this Letter we propose and study a sim- 
ple and general model of evolution of species over a complex 



network. In the model, each species can represent either a 
biotype or a language, a culture or even a consumer product. 
The species compete for space or resources, represented by 
the nodes of the network. For instance, if the species are con- 
sumer products, each node is a potential user and the species 
move from node to node through the network of social rela- 
tionships among users, competing to be adopted by as many 
users as possible. Instead, if we imagine each species as a 
different biotype, the complex network represents the connec- 
tions among spatial environments, and species compete for 
food or energy. The agents of different species move over the 
graph by diffusion and their interaction at the nodes is mod- 
elled by means of competition and replication rules. The aim 
of the present work is to study the combined effects of dif- 
fusion over a complex topology and of competitive selection 
at nodes. Notice that whenever we refer to selection in the 
following, we always intend competitive selection, i.e. com- 
petition for resources among the different species at the same 
node. As we will show in the following, the model, although 
extremely simple, is rich enough to exhibit a large variety of 
patterns over time and a final distribution of the species that is 
intrinsically connected with the structure of the network. 

Let us consider a connected graph with N nodes and K 
links, described by an adjacency matrix A = {a^}, and 
a population of N s different species, labelled by the index 
q = 1, 2, . . . ,N„. Each node of the graph is an environmen- 
tal niche which can host individuals of one or more species, 
and the links represent connections between niches. At ev- 
ery single node of the graph, there is room for each of the 
N s species. We denote with P"(t) the relative abundance 
of the a-th species at node i at time t. Such abundances 
are normalised so that all nodes have the same capacity, i.e. 
£ a Pf(i) = l,Vi,i. We denote with P a (t) = E^PfO) 
the overall abundance of species a in the network, so that 
S Q P a (0 = N . At each time step, the model considers two 
different processes, namely diffusion and competition. In the 
diffusion process, a fraction p (0 ^ p ^ 1) of the individuals 
which are at a given node j move to one of the first neighbours 
of j, let us say i, with a uniform probability: v f J1 . The re- 
maining fraction 1 — p stays at node j. After the diffusion 
process, the selection process takes place, which normalises 
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the number of individuals at each node in order to guarantee 
that ^2 a P"(t) = l,Vi,t. The survival and death of indi- 
viduals at each node is governed by a generalised replicator 
dynamics lfl5i [l6ll . In its simplest version, which takes into 
account an ecosystem with only two species, say X and Y, 
the equations of the replicator dynamics read: 

+ = y(t + l) = 9 ^m y (t). (1) 



where x(t) and y(t) denote the percentages of individuals of 
species X and Y at time i, respectively, while f[x(t)) and 
g[y(t)} are two functions which measure the fitness of each of 
the two species. The interaction between X and Y is ruled by 
the quantity <fr, which plays the role of an environmental limit 
and is fixed to ensure the normalisation x(t) + y(t) = 1 Vi. 
This gives (j> = xf(x) + yg(y), so that (j> = <f>(x,y) is the 
average fitness of the population. The meaning of Eq. ([TJ 
is the following: when the only constraint imposed to the 
species evolution is the environmental limit <f>, their relative 
abundance at the next time step will increase or decrease ac- 
cording to their fitness. As for the fitness functions f(x) and 
g(y) we consider the general case lfl6ll f(x) = bxX 1 ^ 1 and 
g(y) = byy 1 ^ 1 where b x > 0, b y > 0, and the exponent 7 is 
a real number that can be varied to tune the dependence of the 
fitness function of a species on its abundance. The fixed points 
(x* , y*) of Eq. ((TJ, and their stability, depend on the value of 
7. We distinguish three cases: 7 smaller, equal or larger than 
1. For 7 < 1 there are two unstable fixed points (1,0) and 

(0, 1), and one stable fixed point x* = (1 + (bx/by)^) -1 , 
y* = 1 — x*. The two species X and Y will coexist de- 
spite their initial relative abundances and fitness. This case is 
called survival for all. If 7 = 1, there are only two fixed points 
whose stability depends on the respective values of b x and b y . 
When b x > b y then (1,0) is stable while (0, 1) is unstable, 
conversely when b x < b y , (0, 1) is a stable equilibrium while 
(1, 0) in unstable. Independently of the initial distributions of 
the two species, if b x > b y then x will eventually overcome y 
until all individuals of species y are extinct. This case is called 
survival of the fittest. Finally, for 7 > 1 the third fixed point is 
unstable, while (1, 0) and (0, 1) are both stable. In particular 
if the initial condition x(0) is such that x(0) > x*, then x(t) 
will eventually overcome y(t), independently of the value of 
b x and b y , while y(t) will overcome x(t) if x(0) < x* (no- 
tice that x* = 1/2 when b x = by). This super-exponential 
growth always guarantees the survival (and reproduction) of 
the most abundant species, so that the case 7 > 1 is usually 
called survival of the first. In order to implement a strong 
competition among species at a node, in the following we al- 
ways consider 7 > 1. The replicator dynamic can be easily 
extended to N s different species. A super-exponential growth 
is predicted for the A^ -dimensional replicator equations when 
7 > 1, and all the corners of the N s -dimensional simplex are 
stable fixed points. We consider a fitness function of the form 
/(P?) = bf (Pf) T "\ where 7 > 1 and, without any lack of 
generality, we set bf — 1 Vi and Va. In fact, when 7 > 1 only 
the most abundant species will survive, despite the relative 
values of bf (survival of the first). The final model consists of 



the following equations: 
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where i = 1, . . . , N. Eq. © accounts for the diffusion pro- 
cess while Eq. (0 accounts for the selection, with ^ p ^ 1 
and 7 > 1 being the two control parameters of the model. 
The quantity Qf represents the local abundance of species a 
at node i before selection takes place, while P° is the local 
abundance of species a at node i after selection. The dynam- 
ics of the model finally converges to a stationary state with a 
fixed number of surviving species. We have found that this 
number can vary from 1 up to N s , according to the values 
of the two parameters p and 7. Notice that a species a can 
be considered extinct when P a (t) = Pi*(t) < WTtJN' 
where N s (t) is the number of species still present on the net- 
work at time t. This comes directly from the observation that 
when 7 > 1 a species can grow and exponentially reproduce 
on a node if and only if it is the most abundant species on that 
node. If at time t the overall abundance of a species is lower 
than N hs N , then it cannot be the most abundant species at 
any node, and will eventually disappear. A very interesting 
feature of the model is that the number of species at equi- 
librium depends, as expected, on the diffusion (parameter p) 
and on the strength of interaction (the 7 exponent), but it is 
also heavily affected by the topology of the underlying net- 
work. An effective visual representation of the depencence of 
the dynamics on the network structure is the phase diagram 
which reports the number of surviving species as a function 
of p and 7. In fact, the shape of the phase diagram seems 
to be tightly connected with the topological structure of the 
network. We first show how to derive analytically some infor- 
mation on the number of surviving species in the simple case 
of random regular graphs. The fixed points of Eq. (0 and their 
stability can be studied analytically in a mean-field approxi- 
mation. In the mean-field the adjacency matrix of the graph 
is expressed in terms of the probabilty ay = of having 
the edge ay between nodes i and j if fcj and kj are, respec- 
tively, the degree of node i and node j, and K is the number of 
edges in the graph 1I25I1 . Using the mean-field approximation 
for ay in Eq. (f2]i and substituting back in Eq. (0), we obtain 
the following time evolution for the occupation probabilities: 
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We can therefore look for the fixed points: P a (t + 1) = 
P a {t) Va, and check for their linear stability to small per- 
turbations. In the following we consider the case Ns = N, 
i.e. a number of species equal to the number of nodes. It is 
easy to prove that in this case state Si = {P a = l,Va} is 
a stable fixed point for p = 1, and that state S2 = {P a = 
N, P a — Va a} is a stable fixed point Vp. In general, 
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finding all the fixed points of Eqs. ©, for any value of p and 
7, is not an easy task because of the dependencies of the equa- 
tions on the node degrees. A drastic simplification is obtained 
if we make the assumption that all the nodes have the same 
degree k% — k — i.e. when the graph is regular. It is easy 
to verify that for regular graphs in the mean-field approxima- 
tion, the state S3 = {P a = 1, Va; P" = 8i a } is a fixed point 
for all values of p and 7. Notice that in state S3 each node 
contains individuals of only one species, and each species a 
is present only on one node. The edge of the stability region 
for state 5.3 is given by equation: 
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and e is a small perturbation of the fixed point S3. Eq. (|5]l is 
obtained by imposing that S3 is a fixed point, i.e. that relative 
species abundances on the whole graph remain constant over 
time, and then performing a small perturbation on the abun- 
dance of just one species. In particular, we imagine that one of 
the species decreases its abundance by a small amount e, and 
that this amount is uniformly redistributed to the other N — 1 
remaining species, in order to guarantee that the total amount 
of individuals on the network remains constant. Notice that 
Eq.© depends only on p, 7 and N, and does not depend on 
the number of links K, since in the mean-field approximation 
each node is connected to all the other nodes. 




FIG. 1 . Final number of species N a on a random regular graph with 
N = 500 nodes and k — 100 (averaged over 200 realizations), as a 
function of the two parameters p and 7. The yellow area corresponds 
to 500 surviving species, while black area indicates only one surviv- 
ing species. The red line marks the edge of stability for S3 in the 
mean-field approximation. 

The red line drawn in Fig. (Q~|i is the numerical solution of 
Eq. (© for different values of (p, 7) and for e — > on a reg- 
ular network with N = 500 nodes. The state S3 with N 



surviving species is unstable for values of (p, 7) below the 
red line, while it is stable for (p, 7) above the red line. To 
check the validity of the mean-field approximation, we sim- 
ulated the dynamics of Eq. (0 over a regular random graph 
with N = 500 nodes, and a large average degree k = 100, 
for different values of (p, 7), initialising the system in state 
S3. In Fig. (0 we plot as a colour map the number of species 
surviving at equilibrium for different values of p and 7. Yel- 
low regions correspond to 500 surviving species, while black 
regions correspond to only one surviving species. Notice that 
the agreement with the theory is very good: the yellow area 
approximately coincides with the region where S3 is stable, 
while the black zone corresponds to the region where S3 is 
unstable. Notice also that the transition is very sharp and well 
defined in the (p, 7) plane, meaning that the system suddenly 
moves from a state where N s = N to a state where N s = 1. 
To explore the dependence of the dynamics on the degree dis- 
tribution of the network, we have simulated Eq. with initial 
state S3 over three different networks, namely regular lattices 
(RL) with periodic boundary conditions, Erdos-Renyi (ER) 
random graphs, and scale-free (SF) graphs with degree dis- 
tribution P(K) ~ k~ 3 OJ-dl]. All the graphs have been cre- 
ated with the same number of nodes (N = 500) and the same 
number of links ((k) ~ 6). The phase diagrams in Fig. (0 
show the number of species N s remaining at equilibrium on 
the three classes of graphs as a function of the two parameters 
p and 7. We observed that a stationary value of N s is reached 
on these graphs after no more than 200 iterations of Eq. (0, 
and we have checked that this value does not change for at 
least 20000 iterations. Diagrams for ER and SF graphs are ob- 
tained as an average over 200 realizations, even if the fluctua- 
tions from one realization to another are very small. We notice 
that both p and 7 play an important role on the final number of 
species remaining at equilibrium. For small p and large 7 all 
the species survive, with each species remaining in its starting 
node. In fact, when p is close to 0, only a small amount of 
individuals leave their starting nodes and die almost immedi- 
ately after they arrive at neighbouring nodes due to selection. 
This corresponds to the large yellow area (N s = 500) present 
in all the phase diagrams reported in Fig. (0. As the diffusion 
probability p increases, a stronger selection (larger value of 7) 
is needed to prevent the invasion and let all the N s = 500 
species survive. Despite some similarities for small p and 
large 7, the three graphs exhibit different behaviour when dif- 
fusion and selection are such that some species can invade 
neighbouring nodes, and some other species eventually dis- 
appear. According to the values of p and 7, the combination 
of diffusion and competition determines stationary solutions 
with different number N s of surviving species, also with a few 
remaining species, or even just one, as in the black regions. 
The differences between the three graphs are evident from the 
various sizes and shapes of the coloured regions. In both ER 
and scale-free graphs we have cases where only one species 
survives (black regions in Fig. j2(b)\ and in Fig. ( |2(c)| i) and in- 
vades the whole network. However, the black region is much 
larger in the phase diagram of the ER random graphs than 
in that of the scale-free graphs, where a single species over- 
comes all the others only when p ~ 0.95 and 7 ~ 1.2. The 
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FIG. 2. Number of species N s on the graph at equilibrium as a function of the two parameters p and 7. Three classes of graphs with N = 500 
nodes and (k) = 6 were considered: a) regular ring lattice, b) Erdos-Renyi random graph, c) scale-free graph. 



differences between the phase diagrams of ER and SF graphs 
are due to the different degree distribution in the two graphs. 
In fact, the diffusion process naturally favours species start- 
ing at poorly connected nodes, since the average number of 
individuals of such species that will move to first neighbours 
is higher than the average amount of individuals coming from 
highly connected nodes. Hence, species starting at poorly- 
connected nodes have a higher probability to survive and to 
invade neighbours. In the ER graph, these are the few species 
that survive and invade the graph for (p, 7) in the bottom- 
right part of the diagram (purple and black colour), with the 
competition process involving species starting at nodes with 
increasingly large degree, as p decreases. The same consid- 
erations, based on neighbourhood invasion by species starting 
at poorly-connected nodes, hold for scale-free graphs as well, 
with the main difference that in a power-law degree distribu- 
tion the majority of nodes are poorly-connected, while just a 
few hubs have a lot of links. Consequently, species starting 
at hubs will disappear soon, while a large number of species 
tends to survive for a wide range of p and 7. This explains 
why the black region for SF graphs is much smaller than for 
ER graphs, and why at any given point (p, 7) of the phase dia- 
gram, the number of species N s on an ER random graph is al- 
ways equal or lower than the number of species on a scale-free 
graph. The diagram for a regular lattice reported in Fig. ( |2(a)| i 
is very similar to that shown in Fig. (fTJ. These results confirm 
that the degree distribution of the network plays an impor- 
tant role in the extinction and survival of species. We have 
also investigated how the dynamics of the system depends on 
other structural properties of the network, namely the number 
of nodes, the average node degree, the clustering coefficient, 
degree-degree correlations and the presence of communities. 
In the following, we use an alternative method to display the 
information contained in the phase diagrams. Namely, we 
plot the cumulative distribution of the percentage of surviv- 
ing species at equilibrium, over all the couples of parameters 
(p, 7) in the phase diagram. This is useful to compare the 
phase diagrams corresponding to different networks. 



In Fig. ( |3(a)| > we report the cumulative distribution of the 
percentage of surviving species over (p, 7) for three ER ran- 
dom graphs, having average degree (k) = 8 and size N — 
500, N = 1000 and N = 2000, respectively. For the ER 
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FIG. 3. Cumulative distribution of surviving species in ER graphs 
for < p < 1 and 1 < 7 < 6 (averaged over 200 realizations). 
Panel (a) (k) = 8 and N — 500, N = 1000 and N = 2000 nodes. 
Panel (b) N = 1000 and different values of (k). 



random graph with iV = 500 nodes more than 22% of the 
(p, 7) values cause the invasion of the network by one species, 
and N s /N = 1 for more than 45% of the (p, 7) values. 
The percentage of pairs (p, 7) which allow invasion by only 
one species decreases to 20% and to 18%, respectively for 
iV = 1000 and N = 2000, while there is no appreciable dif- 
ference among the three curves for N s /N > 0.4. Therefore, 
for ER random graphs an increase in the network size pro- 
ducees a slight decrease of the area of the phase diagram for 
which we observe invasion by only one species. In Fig. ( |3(b)[ ) 
we compare ER random graphs with N = 1000 nodes and 
different values of (k), namely 6, 8, 10, 15, 20. As the average 
degree increases, the shape of the distribution tends to that of 
a homogenous network. In particular, for (k) = 20, the phase 
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diagram is similar to a stepwise function: only one species 
survives for 25% of the (p, 7) pairs, while all the N s = 1000 
initial species survive for more than 60% of the possible (p, 7) 
values. 

While ER random graphs are characterised by the size of 
the network and by the average degree, networks from the 
real world usually show also non-trivial degree-degree corre- 
lations and a relatively high number of triangles. In Fig. ( |4(a)| i 
we report the cumulative distribution of surviving species at 
equilibrium for the US Airport network ll26ll . This network 
has N = 500 nodes, representing airports, 2980 links, indicat- 
ing flight connections, a degree distribution with a power-law 
tail with an exponent ~ — 1, a clustering coefficient C = 0.62 
and disassortative degree-degree correlations. In the same fig- 
ure we report the cumulative distribution of surviving species 
for a random graph having the same degree sequence of the 
US Airports network. The randomisation washes out all cor- 
relations, so that this network has C = 0.09 and no degree- 
degree correlations. Notice that the US Airports network al- 
lows the survival of a higher percentage of species than its 
random counterpart. More than 85% of (p, 7) values guar- 
antee the survival of more than 50% of the species, and for 
more than 95% of (p, 7) values we observe survival of at least 
40% of the species. Conversely, in the random graph only 
50% of the species survive for 60% of (p, 7) values. This sug- 
gests that even if two networks have exactly the same degree 
distribution, the existence of clustering or degree-degree cor- 
relations favours the survival of a larger number of species. 

We have also found that the presence of modules and 
communities can affect the dynamics of the system. In 
Fig. ( |4(b)| i we show the cumulative distribution of the per- 
centage of surviving species for two standard benchmark net- 
works having a predefined community structure (the Girvan- 
Newman's benchmark (GN) §Zl\\ and the Lancichinetti- 
Fortunato-Radicchi's benchmark (LFR) ||28l0 and for the cor- 
responding random graphs without any community structure. 
All the networks have N = 1000 nodes. The GN benchmark 
is a regular network with k = 12, while the LFR benchmark 
is a scale-free network with P(k) ~ k~ 2 and (k) = 6. The 
two benchmark networks have 4 communities of 250 nodes 
each: on average, 90% of the links of each node are inside its 
community and the remaining 10% of links point to nodes out- 
side the community. The distribution of surviving species in 
the LFR benchmark is very similar to that of a random scale- 
free graph: at least 50% of the species survive for 90% of 
(p, 7) values and more than 70% of the species survive for 
more than 45% of (p, 7) values. Conversely, the GN bench- 
mark has a slightly different behaviour compared to the corre- 
sponding regular random graph. In the regular random graph 
only 1 species eventually invades the network for 22% of pos- 
sible (p, 7) values, while all N s species survive for 65% of 
(p, 7) pairs. For the GN benchmark, instead, more than 20% 
of (p, 7) pairs guarantee the survival of exactly 4 species. 
We have checked that in this case each of the four surviving 
species is confined into one of the four communities, and that 
each node of a given community contains individuals of only 
one species. These results indicate that the existence of com- 



munities in the underlying network can affect the evolution 
of the system, especially in graphs with homogeneous degree 
distributions. 
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FIG. 4. Panel (a): cumulative distribution of surviving species for 
the US Airports network and a random network with the same degree 
distribution. Panel (b): cumulative distribution of surviving species 
for the GN benchmark, the LFR benchmark, a random regular graph 
and a scale-free random graph. The two benchmark networks have 
4 communities of 250 nodes each. 

Summing up, we have found that network structure can 
strongly affect the dynamics of simple diffusion-competition 
processes. A central role in determining the strength of segre- 
gation and the number of surviving species at the equilibrium 
is played by the degree distribution. A network with hetero- 
geneous degree distribution guarantees, for a wider ranges of 
diffusion and selection parameters, the survival of a higher 
number of species compared to the case of a homogeneous 
network. In particular, degree heterogeneity helps to avoid 
the invasion of the network by only one species. We have 
also investigated the effect of other structural properties, such 
as the size of the network, the average degree, the existence 
of degree-degree correlations and community structures. In 
conclusion, the results confirm that the actual structure of the 
network has to be taken seriously into account for the study 
of competitive processes on complex topologies, since small 
differences in the network structure can produce large differ- 
ences in the observed dynamics. Our simple model sheds light 
on the role of the environment in diffusion-competition dy- 
namics, and might find useful to explore how cultures, lan- 
guages, biotypes and competing populations in general may 
survive or get extinct according to the structure of the network 
they live in. 
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